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Abstract 

We theoretically determine the probability distribution function of the net field of the 
random planar structure of dipoles which represent polarized particles. At small surface 
concentrations c of the point dipoles this distribution is expressed in terms of special 
functions. At the surface concentrations of the dipoles as high as 0.6 the dipolar field obey 
the Gaussian law. To obtain the distribution function within transitional region c < 0.6, 
we propose the method based on the cumulant expansion. We calculate the parameters of 
the distributions for some specific configurations of the dipoles. The distribution functions 
of the ordered ensembles of the dipoles at the low and moderate surface concentrations 
have asymmetric shape with respect to distribution medians. The distribution functions 
allow to calculate various physical parameters of two-dimensional interacting nanoparticle 
ensembles. 

The final publication is available at www.epj.org, DOI: 10.1140/epjb/e2012-30057-7 

A considerable interest in studying the electric, magnetic and optical properties of nanocom- 
posites is connected today with two-dimensional systems containing particles. These nanopar- 
ticles may be polarized electrically or magnetically and thus to a first approximation they can 
be regarded as dipoles. Since sufficiently high particle concentrations are widely encountered in 
practice, the long-range dipole-dipole interparticle interaction should be taken into account. In 



particu 
arrays 



ar, the dipolar ferromagnetism was experimentally observed in monolayer nanoparticle 
2). Moreover, the dipole-dipole interactions are substantial for the optical properties 



of nanocomposites 



mi 



The consideration of the dipole-dipole interactions is a complex many-body problem due to 
their long-range and anisotropic behavior. Frequently, this problem is solved using first principle 
numerical simulations. But such a direct computation of the dipole-dipole interactions in an 
ensemble is a resource consuming problem, so we can use the probability distribution function 
of dipolar interaction field in order to allow for interparticle interactions and to derive some 
analytical expressions. After obtaining the probability distribution function, it is possible to 
calculate some physical parameters of a system by means of the methods of the statistical 
mechanics. 

As was shown before, in three-dimensional sample at the low volume concentrations of 
dipoles (spins), the distribution function of the net dipolar field is Cauchy-Lorentzian [sUsI, at 
the large concentrations the distribution becomes Gaussian [7]. The condition for transition to 
the normal law was shown in Ref . ^| . Because of the different geometry, these results cannot be 
applied to the two-dimensional systems. In this study, we determine the probability distribution 
function of the random dipolar field of the flat structure of the spins. 

Let us consider the system of circular particles (nanoislands) which lie randomly in a plane. 
These particles are assumed to have the dipole polarization. We suppose that the particles are 
surrounded by dielectric matrix. At the coordinate origin, net dipolar field projection E onto 
a selected direction falls within range (E, E + dE) with probability 



N 



E 



5> 

fc=i 



dE, 



where £k(Pk, r k) is the field induced by k-th dipole located at position r^, 5 is Dirac's 
delta function. Here, E represents the projection of either electric or magnetic field. Then the 
random field distribution function will be 



W(E)dE 
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N 



Y[r k (p fe ) dp k drdE, 



(1) 



k=l 



where T/~ (p k ) are the distribution functions for dipole moments, dr = r dr dip. After applying 



Markov's method 



lOj , under the assumption that all the r k distributions are identical r k (p* 



r (p)> Pk — P, we obtain characteristic function A(p) = exp [— C(p)\ with 



C( P ) 



J I 1 - ex P (ipQ] r (p) d P dr > 



(2) 



where s is the area of the particle, c = Ns/S is the surface concentration of the dipolar 
polarized particles, S is the area of the sample. It is worth mentioning that C(p) is the 
cumulant-generating function with opposite sign. For the dipolar interaction 

3(r • p)r — pr 2 D 



(3) 



where e m is the dielectric function of surrounding matrix in the case of interacting electric 
dipoles. 

First, we give a consideration to the case of point dipoles which corresponds to limit c — > 0. 
Upon substituting y — 1/r 3 , y\ = l/R 3 , R is the radius of the sample, and neglecting the size 
of the particles, Eq. [2] becomes 

rc (1 - exp(ipDy)) dy ^ f re J 1 - exp(zyipD) ipDT (1/3, -iyipD) 
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where T denotes the incomplete gamma function. For the small values of y\ Eq. |4] takes form 

C(p) = |p| 2/3 ?, 2 '\B 1 + iB 2 signp) - ipE (5) 

with scale parameter 

Bl = 1^(2/3)8 J TlD f /3 dpd ^ (6) 
B2 = -\wFm~s jT\D\ 2/3 signDdpdp (7) 

E = -- s J TDyl /3 dpdp. (8) 

It should be noticed that B 2 and Eq may become zero due to the alternating nature of the 
dipolar fields, for example, when one half of the spins is compensated by another part with the 
opposite direction of the dipoles. On this occasion the system has no ordering. 

The obtained characteristic function corresponds to a stable distribution. It was shown 



skewness 



and median 



in Ref. 



11] that for the infinite two-dimensional ensemble of point sources with some types 



of field dependencies oc l/r M the characteristic function is proportional to exp(— \p\ 2 ^)- Our 



Ising dipoles lying in x-y plane, 
In order to apply inverse Fourier 



result correlate with their findings. For the special case o' 
the similar characteristic function was derived in Ref. [ijjjE 
transform, we split integration over p into two ranges for negative and positive p. Designating 
B = B\ + iB 2 , E = E — Eq, x = B/\E\ 2 ^ 3 we obtain distribution function: 

W{E) = 1 1 sign E [x Bi x 2 + Bi' x 2 } -xAix 2 - Ai' x 2 j ^ exp y - 

1 1 sign E [x* Bi x* 2 + Bi' x* 2 ] + x* Ai x* 2 + Ai' x* 2 j exp + 

( 2 F 2 (1, 1; I, %■ x*) - 2 F 2 ( 1, 1; I, I; x]\ , (9) 



2ttE \ V2 3' 3' J V2 3' 3 



1 An attempt to get the characteristic function in [L.L. Afremov. Dissertation of doctor of sciences, Vladi- 
vostok, Far-Eastern state university (1999)] was unsuccessful due to errors. 



where Ai, Bi, Ai', Bi' are the Airy functions and their derivatives, 2-P2 is the generalized hyper- 
geometric function, the asterisk denotes complex conjugation. In the case of the small values of 
E where the computation of the special functions involves difficulties, the characteristic func- 
tion can be expanded into a series and the integral over p can be taken separately. It gives 
formula: 

TT , /rnN 1 „ (A8B 3 -35E 2 UiE 48£* 3 -35£ 2 UiE } 

= Re (^7^ ^ + —b^ TSF^j ' (10) 

Distribution is asymmetric with respect to its median. For the spin arrangement without 
ordering B2 = and Eq = 0, so function becomes symmetric and simplifies to 

2B ( B 2 B B 2 ^1 IB 3 

W{E) = -jew \ Ai '\EW + \ew m \ew\ exp 3^' 



In this specific case, the distribution was expressed through the Whittaker function in Ref. 1 1 ) 
before: 

W(E) = jL exp ( ^] W i 1 ( ^] . (12) 



Both Eq. [IT] and Eq. [12] are equivalent but the first form is more computationally efficient. In 
the vicinity of zero (E — > 0), Eq. [11] is approximated by formula 

W{E) = 48 "f. (13) 

It should be emphasized that Eq. [9] with alternative parameters B and Eq also describes 
the distribution of interaction fields oc 1/r 4 in the three-dimensional space. 

Fig. [TJ shows functions OH]) and f ill I) in comparison with the Cauchy and Gaussian distribu- 
tions. As well as the Cauchy law, distribution (J9]) has no moments. As one can see from Fig. [1] 
distribution ( II ip is more flat than the Cauchy law. This fact is attributed to the faster reduc- 
tion of ^2 £ with r in two dimensions and the increasing contribution of neighbor spins. The 
main feature of function ([9]) is its asymmetry in contrast to the symmetric three-dimensional 
distribution of the point dipole fields. In fact, the latter arises from the fortunate circumstances 
when C(p) depends only on the integer powers of p. 

Below, we provide the distribution parameters for some spin configurations with identical 
|p| = p. When Na dipoles are oriented perpendicularly to the plane and N/3 dipoles have 
opposite direction (transverse Ising model): 

r( 7 ) = a5(j) + /35(7-tt), (14) 
vr 2 q» 2 / 3 / . r \ cp 2 / 3 (2.023 + 3.504im) 
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Figure 1: Gaussian (1), Cauchy (2), Eq. [TT] (3) and Eq. [9] (4) distributions with a unit area 
under the curves and zero medians. Curve 4 is calculated using B = 0.817 — 0.333z. 

= *£r- (16) 

where m = a — (3 is the order parameter, 7 is the polar angle of p. 

For dipoles lying in the plane and having two possible opposite directions: 

Trcp 2 / 3 R . cp 2 / 3 (1.938 - 1.3692m) 

B \\ = — T^rt ( J i ~ hVStm) « — , (17) 

3V6^/ 3 2s r (2/3) e^ 3 s 

Ji = y |3cosV- l| 2/3 ^ ~ 6.019, 


J 2 = y 1 3 cos 2 ip - 1 1 2/3 sign(3 cos 2 </? - 1) 2.455, 


v ncpm CM 
E ° 11 = "ite> (18) 

When the spins with zero transverse component are arbitrarily oriented in the plane: 

r(7) = L (19) 

2/3 2/3 

B = %R / /"|3cos(2^- 7 ) + cos 7 | 2/3 dip d-f ^ 3.828%-, (20) 



2 5 /337/e e V3 sr (2/3) 







2/3 
Em S 



E = 0. 



These parameters are calculated under the assumption of homogeneous surrounding matrix. 
For other configurations of ambient media they may have more complicated dependence on s m . 
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Further, we investigate the transition to the normal law. At large dipole concentrations c, 
according to the central limit theorem, the stable distribution for the dipolar field switches to 
Gaussian. This can be shown expanding Eq. [2] into series in terms of p up to p 2 . This expansion 
is valid for high c. The projection of net dipolar field onto a selected direction has probability 
distribution function 

W G (E) = ^(-(E^_E )V2^ 

with 



s 



Eo = -- I (r(p)rfpdr (22) 



and variance 



,2 



C 
S 



J er(p)dpdr. (23) 

The integration in Eqs. |22] and [23] over r begins from minimal distance between the particles 
2ro, ro is the radius of the particle. These parameters are readily obtainable for the specific 
configurations of the finite particles. For the Ising model with the transverse alignment of the 
spins: 

Tropin , ncp 2 
E0± ~^ ^-2 S (2r )%^ (24) 

For two possible opposite orientations of the dipoles in the plane: 

TTcpm 2 _ llncp 2 
Eo|| -~2^' a " " 16,(2r )%^' (25) 

and on occasion of the randomly aligned dipoles lying in the plane: 

* = "• ^^k- < 26 > 

It i8 p arti cu. ar ly im po rtan t that the negative sign of E m pe rmlt8 the in-plane dipolar ferro- 
magnetic ordering [13] . 

Making use of the proposed in Ref. 9| comparison of a 2 and the higher moments of spin 
arrangement function r(p), we can evaluate the concentration at which W becomes Gaussian: 

cf£ 3 T(p)dpdr cj i A r(p)dpdr 

3sVNa s ' 12sNa* 

Such an estimate gives that this transition occurs at c > 0.5 for the transverse alignment of the 
spins and at c > 0.6 for the in-plane arrangement of the dipoles. These values are much greater 
than 0.1—0.15 obtained for the three-dimensional sample 9]. Again, it can be accounted for 
the fast saturation of the dipole sums in two dimensions and the relatively large contribution 
of neighbors which differs from the normal law. 
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From the above, we can see that stable distributions (Q and (|2T|) are reached for the two 
extreme cases of the concentrations. Within range of concentrations < c J$ 0.6, it is possible 
to add more terms of expansion ([2]) in powers of p (the cumulant expansion with opposite sign) 

C{p) = E oP + l -a 2 p 2 + ^p 3 p 3 - i/x 4 p 4 + |/i 5 p 5 + • • • , (27) 

where the moments are defined by p n = - J £ n T(p)dpdr and make inverse Fourier transform 
numerically. As well as for Eq and a 2 , it is easy to obtain expressions for the higher moments of 
the simple spin configurations. It is essential that the odd-numbered moments are proportional 
to order parameter m and even-numbered p n do not depend on m due to identity a + = 1. 
The calculations of the distribution functions within the transitional range of the concentrations 
(c = 0.15 and c = 0.3) for in-plane Ising model with m — 1 are depicted in Fig. [2j As well as 
function the negative cumulant expansion gives the asymmetric distribution in the ordered 
state when the odd moments are non-zero. The use of expansion (|27|) up to the 14th degree 
of p shows oscillations of the distribution function which are more pronounced for the lower 
concentrations. For comparison, Fig. [2] illustrates the distributions of the net point dipole 
fields calculated with the same parameters. We must remember that the latter curves are good 
approximation only for c — > 0. The employment of expansion fl27|) benefits when moments 
p n depend on a small parameter which can be applied for further expansion. It is obvious 
from Fig. [2] that distribution differ significantly from the results of the negative cumulant 
expansion at higher c. From this it can be concluded that the model of the point sources is not 
valid at the moderate concentrations. 

It should be underlined that the method based on the negative cumulant expansion is not 
restricted to the two-dimensional geometry and, even, to the dipolar field. This expansion may 
be useful for studying other types of interactions in ensembles with the easy calculation of p n . 

In conclusion, we expressed the distribution function of the random dipolar field in the 
two-dimensional structure in the limit of the low concentrations of dipoles in terms of spe- 
cial functions. This distribution is gently sloping as compared with the Cauchy law obtained 
for three-dimensional sample. At rather high concentrations of polarized nanoparticles, the 
net random dipolar field adhere to the Gaussian law. Within the transitional range of the 
concentrations, the probability distribution function may be calculated after the expansion of 
the negative logarithm of the characteristic function in terms of the spin arrangement func- 
tion moments. The results of this work may be utilized for the calculation of electromagnetic 
parameters of the planar structures with the dipole-dipole interactions. 
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Figure 2: Distributions of the net dipolar field of point dipoles (solid lines), Eq. [91 and finite 
dipoles calculated using expansion (j2T|) in characteristic function up to the 14th degree (dashed 
lines). The geometrical and physical parameters observed in Ref. [2| are used for the simulations. 
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